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We develop a theory for the eigenvalue density of arbitrary non-Hermitian Euclidean matrices. 
Closed equations for the resolvent and the eigenvector correlator are derived. The theory is applied 
to the random Green's matrix relevant to wave propagation in an ensemble of point-like scattering 
centers. This opens a new perspective in the study of wave diffusion, Anderson localization, and 
random lasing. 
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I. INTRODUCTION 

Random matrix theory is a powerful tool of modern 
theoretical physics [l]. First introduced by Wishart [2 
and then used by Wigner to describe the statistics of en- 
ergy levéis in complex nuclei random matrices are 
nowadays omnipresent in physics [4]-[7]. The majority 
of works — including the seminal papers by Wigner [3] 
and Dyson [8 — deal with Hermitian matrices. Her- 
mitian matrices are, of course, of special importance in 
physics because of the Hermiticity of operators associ- 
ated with observables in quantum mechanics. However, 
non-Hermitian random matrices also attracted consider- 
able attention in particular because they can be 
used as models for dissipative or open physical systems 

[EHH]. 

A special class of random matrices are the so-called 
Euclidean random matrices (ERMs) [15 . The elements 
Aij oí di N X N Euclidean random matrix A are given by 
a deterministic function / of positions of pairs of points 
that are randomly distributed in a finite región V of Eu- 
clidean space: Aij = /(r^,rj), i = 1,...,A/'. Hermitian 
ERM models play an important role in the theoretical de- 
scription of supercooled liquids p!5HT8) , disordered super- 
conductors relaxation in glasses and scalar phonon 
localization [20] . They have been used as a playground to 
study Anderson localization [21]. A number of analytic 
approaches were developed to deal with Hermitian ERMs 
[15-22^ . Non-Hermitian ERMs appear in such important 
physical problems as Anderson localization of light [23] 
and matter waves [24 , random lasing [25 , propagation 
of light in nonlinear disordered media [26 , and collective 
spontaneous emission of atomic systems [271 El] • How- 
ever, no analytic theory is available to deal with non- 
Hermitian ERMs and our knowledge about their statisti- 
cal properties is based exclusively on large-scale numer- 
ical simulations [22r>26j . The principal difñculties that 
one encounters when trying to develop a theory of non- 
Hermitian ERMs stem from the nontrivial statistics of 
their elements and the correlations between them. Both 
are not known analytically and are often difñcult to cal- 
cúlate. This is in contrast with the works [12-14 where 
the joint probability distribution of the elements of the 
random matrix under study is the starting point of anal- 
ysis. 



In the present paper we develop an analytic theory for 
the density of eigenvalues of an arbitrary non-Hermitian 
ERM in the limit of large matrix size (TV oo). Par- 
ticularly simple results are obtained for the borderline of 
the support of eigenvalue density on the complex plañe. 
We illustrate the power of our approach by applying it to 
the 'random Green's matrix' — a matrix with elements 
given by the Green's function of the scalar Helmholtz 
equation — that previously appeared in Refs. ¡22H28] but 
was studied only numerically up to now. We discuss the 
link that exists between our calculation and the theory of 
wave scattering in disordered media as well as the local- 
ization properties of eigenvectors of the random Green's 
matrix. 



II. FOUNDATIONS OF THE NON-HERMITIAN 
RANDOM MATRIX THEORY 

The density p{A) of eigenvalues A of any random NxN 
matrix A can be obtained from the resolvent 

If A is Hermitian and A are real, one conveniently ex- 
pands g{z) in series in 1/z in the vicinity of \z\ oo, 
performs the calculation using diagrammatic or any other 
approach, and use the result obtained after the resum- 
mation of the series at all z to obtain 

p{A) = - - lim Im^(A + ze). (2) 

TT e^0+ 

For a non-Hermitian matrix A, however, A are complex 
and g{z) loses its analyticity inside a two-dimensional 
domain V on the complex plañe where A are concen- 
trated. Thus, g{z) íot z e V cannot be assessed by ana- 
lytic continuation of its series expansión in the vicinity of 
\z\ oo. A way to circumvent this problem is to double 
the size of the matrix and to work with a new 2N x 2A^ 
matrix 

A^-{i¡,) (3) 

for which the generalized resolvent matrix 



2 



is safely equal to its series expansión [10]. Here Ttn 
denotes the block trace of a 2N x 2N matrix [see Eq. 
( A2 ) of Appendix [A] for the definition] and 
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(5) 



The resolvent g{z) can be found from the diagonal ele- 
ments of the 2x2 matrix obtained by taking the hmit 
e ^ 0+ in Eq. ^ : 



hm G(Z,) = 

6^0 + 



g{z) c{z) 
c{z) g{zY 



(6) 



and the density of eigenvalues A inside its support V on 
the complex plañe is [10 



TT dz* 



(7) 



with the standard notation d/dz* = ^{d/dx -\- id/dy) 
ÍOT z = X -\- iy. The off-diagonal elements of G yield the 
correlator of right \Rn) and left eigenvectors of A 



N 
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Nc{zf 



(8) 



III. NON-HERMITIAN EUCLIDEAN RANDOM 
MATRIX THEORY 

All above applies to any non-Hermitian matrix A. Let 
US now make use of the fact that A is an ERM with el- 
ements Aij = /{vi^Vj) = (r¿|Á|rj). Here the N points 
Ti are randomly distributed inside some región V of d- 
dimensional space with a uniform density p = N/V^ 
and we introduced an operator A associated with the 
matrix A. A useful trick consists in changing the ba- 
sis from {vi} to {^|Ja} which is orthonormal in V [22 . 
In a rectangular box, for example, li/ja) = \^a) with 
(r|kc^) = exp(ikc^r)/V^ [22 . For arbitrary V we have 



a = hth\ 



(9) 



where Hio, = (r¿|V^c.)/v^ ^c./3 = p {^l^a\A\^lJ^) [34 . 
The advantage of this representation lies in the separa- 
tion of two different sources of complexity: the matrix H 
is random but independent of the function /, whereas the 
matrix T depends on / but is not random. Furthermore, 
if we assume that (Hia) = O, which in a box is obeyed for 
all a except when k^^ =0, we readily find that Hi^ are 
identically distributed random variables with zero mean 
and variance equal to 1/A^. We will assume, in addi- 
tion, that Hia are independent Gaussian random vari- 
ables. This assumption largely simplifies calculations but 
may limit applicability of our results at high densities of 



points p, at least for certain types of Euclidean matrices, 
as we will see later. 

In Appendix |Aj using the diagrammatic expansión of 
the self-energy matrix S(Ze) = — G{Ze)~^^ we show 
that, due to the representation ([9| and the Gaussian 
statistics of ií, in the limit of large ^(^e) involves 
only planar rainbow-like diagrams [10 . Summation of 
these diagrams yields coupled equations for operators 
¿11 and ¿12 that give the elements Su = TtÍ^h/N and 
Si2 = TrÉi2/A^ of the 2 X 2 matrix S = lim,^o+ ^(^e): 



Éii = (l+^Éii+cÉi2)T, 
Éi2 = (cÉii+^*Éi2)Tt, 



(10) 
(11) 



where T = pÁ. After some algebra, these equations lead 
to two self-consistent equations for the resolvent g{z) and 
the eigenvector correlator c{z): 



-Tr 



(1 -í/*ft)f 



N (l-^*Tt)(l-í/T)-c2TtT 



-Tr- 



c2 TV (1 -^*Tt)(l -^T) -c^TtT 



,(12) 



(13) 



Because c{z) should vanish on the boundary 5V of the 
support of the eigenvalue density P, equations for z ^ 5V 
follow: 



z=-^ ¿Tr^, 
g N 



(14) 
(15) 



where 5' = T/(l - 

Equations ([12|, (13), (14) and (15) are our main re- 
sults. An equation for the borderline of the support of 
the eigenvalue density of a non-Hermitian ERM A on 
the complex plañe z = K follows from Eqs. (14) and (15) 



gf). 




upon elimination of g. The density of eigenvalues A in- 
side its support V can be found by solving Eqs. 
(13) with respect to g{z) and then applying Eq. 
analysis includes the result for Hermitian ERMs as a spe- 
cial case: if A is Hermitian, then S is diagonal and the 
support of the eigenvalue density shrinks to a segment on 
the real axis. Equation ([T4| then allows one to solve for 
g{z). This result for Hermitian matrices coincides with 
the one found in Ref. [22 using a different approach. 
The solution of Eqs. 



12L (13), (14) and (15) for a 
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given matrix A is greatly facilitated 
of the basis in which traces appearing in these equa- 
tions are expressed. In addition to {r} and {k^^}, a bi- 
orthogonal basis of right \Ra) and left \Ccx) eigenvectors 
of T can be quite convenient. The right eigenvector \lZo,) 
obeys 

(r|f |7^«) =p j dV/(r,r07^«(rO = /ic.7^c.(r), (16) 
Jv 

where /íq, is the eigenvalue corresponding to the eigen- 
vector \lZcy)' The traces appearing in Eqs. (14) and (15) 



3 



can be expressed as 

Tr^ = J2{c^\s\na) = J2Y 



5Ma 



TrSS'' = J2 



a,/3 



ll^ll*p{C^\Cp){Tlp\Ra) 

(1- giioc)ií- gi^lí)* 



(17) 
(18) 



respectively. 



IV. EXAMPLE: RANDOM GREEN'S MATRIX 

Let US now illustrate the power of the above general 
analysis on the example of the random Green's matrix 



(1 



Sij) 



ex.p{iko\Ti - Tj\) 

ko\Ti-Tj\ 



(19) 



where ko = 27r/Ao and Aq is the wavelength. We as- 
sume that the N points are chosen randomly inside 
a three-dimensional {d = 3) sphere of radius R. This 
non-Hermitian ERM is of special importance in the con- 
text of wave propagation in disordered media because 
its elements are proportional to the Green's function of 
Helmholtz equation, with r¿ that may be thought of as 
positions of point-hke scattering centers. It previously 
appeared in Refs. [22-28J, but was studied only by ex- 
tensive numerical simulations, except in Ref. [28] where 
analytic results were obtained in the infinite density hmit. 

For each reahzation of the random matrix (19), its 
eigenvalues obey [22] 



N 

E 

n=l 



An = O, ImAn > 



(20) 



Very generahy, the eigenvalue density of the matrix de- 
fined by Eq. (19) depends on two dimensionless param- 
eters: the number of points per wavelength cubed pAg 
and the second moment of |A| calculated in the limit of 
low density: (|Ap) = 7 = 9N/8{koR)'^. Even though the 
latter result for (|Ap) can be rigourously justified only in 
the limit of low density pX^ <C 1, it holds approximately 
up to densities as high as pAg ~ 100. We will see from 
the following that the two parameters pAg and 7 control 
different properties of the eigenvalue density. 



A. Borderline of the eigenvalue domain 

We first focus on the borderline of the support of eigen- 
values which is easier to visualize. In Fig. [l]we present 
a comparison of the solutions of Eqs. ([T4|) and (ITsj) with 



results of numerical diagonalization of the matrix (19) 



for koR 1. At low density pX^ < 10, a sufíiciently 
accurate solution of Eqs. ( [m] ) and ( [Ts] ) can be obtained 
in the |r) representation, in which 



{r\S\r')c^p 



exp(ÍK|r — r'l) 
fco|r-r'| ' 



(21) 



with hv{g) = /cq vT^r^pA^T^TT^- Appendix [b| we show 
that this leads to a borderline equation 



where 



h{x) 



|Ap = 27/1 [2lmi^{l/A) R], 



3 - + g^3 _ ^(^i ^ 2x)e- 



-2x 



For pAg ^ 10, a simpler equation 



|A|^ 2-fh -87 



ImA 



(22) 



(23) 



(24) 



yields satisfactory results as well. For 7 <C 1, the den- 
sity of eigenvalues is roughly uniform within a circular 
domain of radius see Fig. [ija). The domain grows 

in size and shifts up upon increasing 7. At 7 > 1 it starts 
to 'feel' the 'wall' ImA = — 1 and deforms [Fig. [TJb)]. 

The approximate equation ( [22] ) for the borderline of 
the support of eigenvalue density yields a closed line on 
the complex plañe until pAg — 30, after which the line 
opens from below. This signáis that an important change 
in behavior might be expected at this density. And in- 
deed, we observe that a 'hole' opens in the eigenvalue 
density for pAg ^30. As we see in Fig. [l]^c), this hole is 
perfectly described by our Eqs. (14) and ( [l5| ) which we 
now solve in the basis of eigenvectors of the operator T. 
Eigenvalues and eigenvectors of T can be found analyti- 
cally [28]. As we discuss in Appendix[Cj this allows for an 
exact solution of Eqs. (ImI) and (15). Finally, at very high 



density the crown formed by the eigenvalues blows up in 
spots centered around the eigenvalues fia of T, as we 
show in Fig.[TJd). When the density is further increased, 
the eigenvalues A^ of A become equal to the eigenvalues 
Pq o f T. They then fall on the circular line given by Eq. 
(C7) and the problem looses its statistical nature. As 



follows from our analysis, the parameter 7 controls the 
overall extent of the support of eigenvalue density V on 
the complex plañe, whereas its structure depends also on 
the density pAg. At fixed 7, V goes through a transition 
from a disk-like to an annulus-like shape, and eventu- 
ally splits into múltiple disconnected spots upon increas- 
ing pAg. The transition from disk-like to the annulus- 
like shape is reminiscent of the disk-annulus transition in 
the eigenvalue distribution of rotationally invariant non- 
Hermitian random matrix ensembles ¡TI] . 

An important additional feature of the numerical re- 
sults in Fig. [1] that is not described by our Eqs. (14) 
and (15) is the eigenvalues that concéntrate around the 
two hyperbolic spirals, |A| = 1/ argA and its refiection 
through the origin. These spirals correspond to the two 
eigenvalues ±^12 of the matrix ([l9| for A/^ = 2 \22í\23\. 
The eigenvectors corresponding to these eigenvalues are 
localized on pairs of very cióse points. From numeri- 
cal results for N < 10^, we estímate their statistical 
weight to be important at large densities, of the order of 
1 — const/(pAo)^ with p ^ 1. This is consistent with the 
estimation of the number of subradiant states in a large 
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FIG. 1: Density plots of the logarithm of eigenvalue density of the N x N random Green's matrix (19) obtained by numerical 



diagonalization of 10 realizations of the matrix for N = 10^. The sohd hnes represent the borderhnes of the support of eigenvalue 
density following from Eq. (22) in panels (a) and (b) and from Eqs. (C12) and (C13) of Appendix[C]in panels (c) and (d). The 
dashed hnes show the diffusion approximation ( 34 ) . 



atomic cloud by Ernst [2T. At large densities, the ab- 
solute majority of the lacking eigenvalues fall very cióse 
to the axis ImA = —1, in the 'gap' that opens in the 
eigenvalue distribution following from our theory on the 
left from ReA = O [see Figs. [TJc) and (d)]. The lack of 
the spiral branches of p{A) in our theory can be traced 
back to the assumption of statistical independence of el- 
ements of the matrix H in Eq. ([9|. It does not affect the 
excellent agreement of the borderline of the rest of the 
eigenvalue domain with numerical results. 



resonant point-like scatterers. The latter problem is de- 
scribed by the Helmholtz equation associated with a fic- 
titious Hamiltonian 



N 



(25) 



i=l 



The retarded free-space Green's function corresponding 
to v{ko) = O, 



B. Mapping to the scattering theory 

We now want to introduce an interesting mapping be- 
tween our results for the random Green's matrix (19) 



and the problem of múltiple scattering of waves by N 



1 



^2 + + V2 ' 



is simply proportional to the matrix (19): 



{go)ij = {ri\úo\rj 



ko_ 
47r^ 



(26) 



(27) 
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Expanding the Green's function 

5= ' 



kl+ie- H 



in Born series, we get 
Q 



(28) 



(29) 



where t is the scattering matrix of an individual scatterer 
defined by ¡SO] 



tó(3)(f-r,) 



v{ko) + v{ko)5^''\v - r¿)áoí] S^^'Kv - r¿). 

(30) 



At , the intensity of a wave emitted by a point source 
located at Yj is lij = \Gij\^^ where Qij = (r¿|í?|rj). Let 
US introduce I{t) = ^í-lj ^iji where we emphasize that / 
depends ont. It can be readily written as 



m = Tr 



(í_g-i)(í_g-i)t 



(31) 



This is to be compared with the expression for the corre- 
lator of right and left eigenvectors of an arbitrary matrix 
A fohowing from Eq. (|6|: 



= — hm — Tr ( — — r - 

^ ^ e^o+ N \{z- A){z - A)t + e2 

For A = Qq^ and z = t we thus have 



c(í) = 



(32) 



(33) 



This should become different from zero when t enters the 
support of the eigenvalue density of Gq^ or, equivalently, 
when 1/t enters the support of the eigenvalue density of 
Go' The only way to obtain c{t) 7^ O for e ^ 0+ is to 
make I{t) diverge. In the framework of our linear model 
of scattering, this can be achieved by realizing a random 
láser |31j. We thus come to the surprising conclusión 
that finding the borderline of the support of the eigen- 
value density p{A) oí the N x N Green's matrix (19) is 
mathematically equivalent to calculating the threshold 
for random lasing in an ensemble of N identical point- 
like scatterers with scattering matrix t = —Air/koA. In 
the diffusion approximation, for example, the threshold 
of such a random láser can be found as in Ref . [32] . This 
leads to the following equation for the borderline [35] : 



IM 



87 



Vi + ImA 1 



lAI- 



|A|2+47 



(34) 



We show this equation in Figs. [TJa) and (b) by dashed 
lines. As expected, it gives satisfactory results only in 
the weak scattering regime pAg ^ 10 and at large optical 
thickness b = 2R/Í = 167/3|Ap > 1, where _£ = At: / pW 



is the mean free path. In contrast, our Eqs. ( 14 ) and ( 15 ) 
apply at any pAg and h. These equations can therefore 
serve as a benchmark for theories of múltiple scattering. 
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FIG. 2: Cuts of the eigenvalue density p(A) of the Green's 
matrix (19) along the imaginary axis Re A = 0. Numerical 



simulations (symbols) are compared with our analytical solu- 
tion (35). 



C. Eigenvalue density 

Let US now analyze the shape of the eigenvalue density 
p(A) inside its support V using Eqs. (12) and (13). Very 



generally, p(A) is roughly symmetric with respect to the 
line Re A = O and decays with ImA. A particular feature 
of p(A) that was studied previously is the behavior of the 
marginal probability density of ImA. Pinheiro et al. ¡23] 
observed p(ImA) oc l/(ImA + l) in numerical simulations 
at high density and conjectured that it was a signature 
of Anderson localization of waves in the corresponding 
point-scatterer model. To test this conjecture, we analyze 
p(A) at low densities pAg < 1, for which no Anderson 
localization is expected. An approximate solution of Eqs. 
( 12 ) and ( 13 ) in this regime can be obtained by neglecting 



the term c^f^f in their denominators: 



— TrS^ 



(35) 



Traces in this equation can be explicitly calculated us- 
ing Eq . ([21] ) valid at low densities, as we show in Ap- 
pendix |B[ The eigenvalue density p{A) is then found by 
applying Eq. ([t]). In Fig. [2] we show cuts of p{A) along 
the imaginary axis Re A = 0. We clearly observe that 
p(ReA = O, ImA) decays as l/(ImA + 1), even though 
the density of points pAg is too low to bring the system 
to the Anderson localization transition. However, the 
power-law decay becomes clearly visible in the marginal 
distribution p(ImA) only when the support of p(ImA) is 
sufíiciently wide, i.e. for 7 > 1. Otherwise, it is 'spoiled' 
by the circular shape of the support of p{A) and p(ImA) 
follows the Marchenko-Pastur law [22]. Because the con- 
dition 7^1 can be obeyed at any, even very low density 
by just increasing the number of points it seems that 



no direct link can be established between the power-law 
decay of p(lmA) and Anderson locahzation. 



D. Anderson localization 

It should be stressed here that Anderson locahzation 
— the localization of eigenvectors in space due to disor- 
der — is a property of eigenvectors \Rn) of the matrix 



(19), whereas our study in this paper concerns its eigen- 



values An. It is not clear a priori if any sign of Anderson 
localization should (and could) be visible in the density 
of eigenvalues p{A). To elabórate on this issue, we an- 
alyze the eigenvectors of the matrix (19). To determine 
if an eigenvector \Rn) is localized or not, we compute its 
inverse participation ratio (IPR): 



IPRn 



2 • 



(36) 



An eigenvector extended over all N points is character- 
ized by IPR ~ l/A/", whereas an eigenvector localized on a 
single point has IPR = 1. The average valué of IPR corre- 
sponding to eigenvectors with eigenvalues in the vicinity 
of A can be defined as 



IPR(A) = -i- / ^ IPR„ 5\A - A„) 



N 



p{A) 



(37) 



\n=l 



where averaging is over all possible configurations of TV 
points in a sphere. Our numerical analysis of the aver- 
age IPR defined by this equation reveáis the following 
scenario. At low density pX^ < 10, IPR ^ 2/N for all 
eigenvectors except those corresponding to the eigenval- 
ues that belong to spiral branches in Fig. [TJa) and (b) 




FIG. 3: Density plot of the logarithm of the average inverse 
participation ratio of eigenvectors of the Green's matrix (19). 



To obtain this plot, we found eigenvalues of 10 different ran- 
dom realizations of 10^ x 10^ Green's matrix numerically, com- 
puted their IPRs using Eq. (36 ), and then determined IPR(A) 
by integrating Eq. (37) over a small área (AA)^ around A, for 
a grid of A's on the complex plañe. 



the limit of ^ oo. These equations allow us to analyze 
the borderline of the support of eigenvalues A by looking 
for a contour on the complex plañe on which c{z) = O, 
as well as the full probability density p{A) inside this 
contour by solving for g{z). To give an example of ap- 
plication of our general results to a particular physical 
problem, we studied the eigenvalue density of the ran- 



dom Green's matrix (19). An entry Aij oí this matrix 



is equal to the Green'smnction of the scalar Helmholtz 



for which IPR \. These states are localized on pairs equation between two points r¿ and Vj chosen among N 

points randomly distributed in a sphere. We showed that 
finding the borderline of the support of the eigenvalue 
density of the Green's matrix is mathematically equiva- 
lent to calculating the threshold for random lasing in an 
ensemble of N identical point-like scatterers. Finally, we 
discussed manifestations of Anderson localization in the 
properties of this matrix and challenged the link that was 
previously proposed between Anderson localization and 
the power-law decay of the marginal probability density. 



of points that are very cióse together and correspond to 
proximity resonances [23 that do not require a large op- 
tical thickness to build up. The prefactor 2 in the result 
for IPR of extended eigenvectors is due to the Gaussian 
statistics of eigenvectors at low densities. For pAg ^ 10, 
IPR starts to grow in a roughly circular domain in the 
vicinity of A = O and reaches máximum valúes ~ 0.1 at 
pAq 30 [see Fig. [s]. Contrary to common belief [23 , 
neither localized states necessarily have ImA cióse to — 1, 
ñor states with ImA — 1 are always localized, as can be 
seen from Fig. [sj For pAg > 30, the localized states start 
to disappear and a hole opens in the eigenvalue density. 
It is quite remarkable that the opening of the hole in 
p{A) proceeds by disappearance of localized states (i.e., 
of states with IPR > l/N). 
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V. CONCLUSION 

We derived equations for the resolvent g{z) and the 
correlator c{z) of right and left eigenvectors of an arbi- 
trary N x N non-Hermitian Euclidean random matrix in 



Appendix A: Derivation of self-consistent equations 
for the resolvent and the eigenvector correlator 



The purpose of this Appendix is to derive Eqs. (12) 



and (13) of the main text. We start by expanding the 
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2x2 resolvent matrix G{Z^) defined by Eq. Q in series 
in 1/Z, = {1/Z,)^1n: 



1 

TV 



Tr 



N 



Gi2 

1 



1 



(Al) 



where the averaging (...) is performed over the ensemble 
of matrices H entering the representation ([9| of the ma- 
trix A. The block trace TtnX of an arbitrary 2N x 2N 
matrix X is defined by separating X in four N xN blocks 
Xii, Xi2, X21, X22 and taking the trace of each of the 
latter separately: 



TtnX 



Tr 



Xii X12 
X21 X22 

[TiXn TrXi2 
V TrX2i TrX22 



(A2) 



As explained in the main text, we assume that H has 
independent identicahy distributed complex entries that 
obey circular Gaussian distribution. Using the properties 
of G auss ian random variables, the result of averaging in 
Eq. (Al) can be expressed through pairwise contractions 



1 

TV 



i/3 



(A3) 



To evalúate efficiently the weight of different terms that 
arise in the calculation, it is convenient to introduce dia- 
grammatic notations. First, the matrices ií, ií^, A and 
A^ will be represented as shown in Fig. [ij 



H = J 



H^= i L 



V\.fJL 



O 



W. f\¡L 



FIG. 4: Diagrammatic representations of the matrices H, 
H\ A = HTH^ and A^ . FuU and dashed Unes propágate in 
the bases {r^} and {V^a}, respectively. 



The 'propagator' 1/ will be depicted by 



1 



-\z\^ 



1112 
2 r 2 2 



(A4) 



Each contraction ( |A3[ ) brings a factor and each 

loop corresponding to taking the trace of a matrix brings 
a factor see Fig. [5] 

In the limit N ^ 00^ only the diagrams that contain 
as many loops as contractions will survive. These dia- 
grams are those where full and dashed lines do not cross. 
Therefore, the leading order expansión of the resolvent 
(Al) involves only diagrams which are planar and look 
like rainbows. Such diagrams appear, for example, in 



íí^'\ = ír^\ = i/N 

= Nx, í[.X^ = TrX 

FIG. 5: Diagrammatic notation for pairwise contractions 
(A3) and loop diagrams for any scalar x in the basis {fí}, 
and for any operator X in an arbitrary basis {ipa}- 




FIG. 6: Diagrammatic expansión of the two independent 
elements of the matrix G(Ze)- 



Fig. [6j where we show the beginning of the expansión of 
the two independent elements of G{Z^). 

In the standard way, rather than summing up the dia- 
grams for the resolvent, we introduce the 2x2 self-energy 
matrix 



S(Z,) = Z,-G(Z,)-i 



_ f Sfi S 

^12 ^11 



12 



(A5) 



It is equal to the sum of all one-particle irreducible dia- 
grams contained in 



1 

iV' 



Tryv < A 



■ (A6) 



The first dominant terms that appear in the expansión 
of the two matrix elements Ya\y and represented 
in Fig. I?! In the two series of Fig. [7|we recognize, under 
a pairwise contraction, the matrix elements G\^ and G\2 
depicted in Fig. |6| as well as the two operators and 
defined in Fig. [sj 
Equations obeyed by the operators ¿11= lime^o+ ¿11 
and £12 = lim£^o+ ¿12 obtained after summation of 
all planar rainbow diagrams in the expansión of Fig. [7| 
and taking the limit e ^ 0+. The diagrammatic repre- 
sentation of these equations is shown in Fig. [Üj Equations 
([To| and (11) of the main text follow after application of 
'Feynman' rules defined in Fig. [5] Furthermore, as fol- 
lows from Eq. (|6]) and the definition of the self-energy 
matrix, in the limit e ^ 0+, ^ and c are simply related 
to En = TrÉii/TV and S12 = TrÉi2/TV by 



-S12 ^* 



-^12 

-Sil 



(A7) 



Elimination of the self-energy S from Eqs. ( 10 ), ( 11 ) and 



(A7) yields Eqs. (12) and (13) of the main text. 




FIG. 7: Diagrammatic expansión of the two independen! elements of the self-energy Braces with arrows denote parts 

of diagrams that are beginning of diagrammatic expansions of the quantities which the arrows point to. 



ye — (¡ — \i 'A 

FIG. 8: The elements Sii and of the matrix S(^e) can 
be written as traces of operators ¿ii and £12 that appear in 
Fig-I?! Sil = Tráíi/AT and = TÁl^jN . 



^12 ~~ 



FIG. 9: Coupled equations for the operators Su and S12 
that define the self-energy S = lim£^o+ ^{Ze)- Here g = 
limg^Q+ Gil and c = limg^Q+ GI2 [see Eq. mí]. 



Appendix B: Approximate solutions for the 
borderline of the eigenvalue domain and the 
eigenvalue density at low density 

Let US show how an explicit equation for the border- 
line of the support of eigenvalue density of the random 
Green's matrix (19) — Eq. — can be derived in the 



low-density limit. On the one hand, traces appearing in 
Eqs. (14) and (15) in the |r) representation read 



TrS' = Tr 



Tr 



(f + gfS 



g /^^d^rdVT(r,rO^(r^r), (Bl) 



/X 



d^r d^r'|5(r,rO 



(B2) 



where T(r, rQ = p{r\Á\r') = pexp{iko\r - r'|)/A:o|r - r'| 
and in Eq. (|bi|) we used the fact that TrT = pTri = O, 



as follows from Eq. ([19]). On the other hand, 6'(r, r') = 
(rl^lr^) obeys 

5(r, rO = T(r, r^) ^ g í d^r''T(r, r^^)S{r^\ r^, (B3) 
Jv 



as follows from the definition of S. Noting that 

4:7Tp 



(Ar 



k¡ + ie)T{r,r') 



ko 



(B4) 



where e ^ 0+, we apply the operator Ar + /cq + ie to Eq. 
(B3) and obtain 



Ar5(r,r') + fc2 



l+g0nv{v) + ie 



S{r,r') 



-^á(3)(r-r'), (B5) 



where ny(r) = 1 for r G F and O elsewhere. In the 
limit of low density pX^ ^ O, an approximate solu- 
tion of this equation is obtained by neglecting 'reflec- 
tions' of the 'wave' 6'(r, r') on the boundaries of the 
volume V and thus setting ny(r) = 1 everywhere. 
This yields 6'(r, r') :^ pexp(zAí|r — r'\)/ko\r — r'\ with 
K{g) = koVl^gpXl/2^^. 



In order to evalúate the integráis (Bl) and (B2), we 



will make use of the following auxiliary result: 
rr d^r c\^r' 

// ir-í^/(l^-^'l)=24 / dxf{2Rx)s{x)x\ 

J Jv{R) V V Jq 

(B6) 

where / is an arbitrary function, V{R) = 47ri^^/3, and 
s{x) = 1 — 3x/2+x^/2. To derive this equation, we define 
new variables x = (r — y')/2R and y = (r + v')/2R. The 



9 



conditions r < R^r' < R become + y 
with O < t < 1, so that 



2xyt < 1, 



// 

J JV{R) 



V V 



-(•••) 



- / d^x / di/ dy-2 
^ ^y(i) Jo Jo 



yX-), 



where yM{t^x) = + (í^ — l)x 



(B7) 

-tx. Evaluation of all 



integráis except one in Eq. (B7) leads to Eq. (B6). 



We now plug t he ex plicit exp ressions for T (r, r ^) and 
6'(r,rO into Eqs. ^ and (|B2| and use Eq. (|B6|. This 
yields 



TrS' = 2-fNgh[-m{g)R - ik^R], 
TrSS^ = 2-iNh[2lmK.{g)R], 



(B8) 
(B9) 



with 



h{x) 



si ^us{u)e 



-2ux 



jl áus{u) 



1 



- , [3-6x^ + -3(l + 2x)e-^^], (BIO) 

and 7 = 9N/S{koR)'^ . In the low density hmit, the latter 
is equal to the second moment of the absolute valué of 
A: 7 = (|Ap). We checked numerically that even at 
higher densities (at least, up to pAg ~ 100), 7 is still a 
good approximation for (|Ap) and henee a meaningful 
parameter. 

In the low-density limit, g can be eliminated from Eqs. 
( [l4| and (15) by neglecting TtS/N in Eq. (14) and sub- 
stituting g = l/z into Eq. (|B9|. This yields Eq. ([22]) and 
then Eq. ( [24| ), if the argument of the function h in Eq . 
(22) is expanded in series in pX^. By comparing Eq. ( |24| ) 
with the exact solution obtained in Appendix[Cj we con- 
clude that it is valid up t o de nsities as high as pX^ o:^ 10. 

Finally, Eqs. (|b8| and (|b9| for Tr6' and TtSS"^ can be 
used to find the resol ve nt g{z) using Eq. (35) and then 
the density of eigenvalues p(A) using Eq. 



Appendix C: Exact solution for the borderline of the 
eigenvalue domain at any density 



In this Appendix we show how Eqs. (14) and (15) can 
be solved exactly using the bi-orthogonal basis of right 
\Tía) and left \Ca) eigenvectors of T. These eigenvectors 
obey T\1Za) = Pa\^a) and T^C^) = f^al^a)- this 
basis, Eqs. (14) and (ÍIS]) read 



1^1 



(1 -^Ma)(l 



1 _ 1 ^ /J^aí 



(C2) 



where, similarly to the derivation in Appendix [B| we 
made use of the fact that TrT = O and therefore TiS = 
gTiTS. The problem essentially reduces to solving the 
eigenvalue equation 



p / d'^r 
Jv 



3 ,exp(ifco|r 



A:o|r - 



-7^,(r0 = /i«7^«(r), (C3) 



where r e V. As follows from Eq. (B4), 7^a(r) is also 
an eigenvector of the Laplacian operator, Ar7^Q;(r) = 
— /T:^7^Q,(r), with k.^ = ^(l/^a)- Iii a sphere of radius 
i?, using the decomposition of the kernel of Eq. ( C3 ) in 



spherical harmonics, it is quite easy to find that [28 

7^c.(r) = Uimpir) = Aipji{^ipr)Yi^{0, (/>), (C4) 

where O and (j) are the polar and azimuthal angles of the 
vector r, respectively, ji are spherical Bessel functions of 
the first kind, Yi^ are spherical harmonics, Aip are nor- 
malization coefficients, and a = {/,m,p}. Furthermore, 
coefíicients tzip obey [28^ 



ko 



ji-i{i^ipR) hl^\koR) ' 



(C5) 



where /¿[^^ are spherical Hankel functions. Integer p la- 
bels the different solutions of this equation for a given /. 
Henee, eigenvalues pip = pAo/27r^(/í^^//co — 1) are (2/ + 1)- 
times degenerate (m G [—l^l])- 

In the limit koR 1, we can use as ymp totic expres- 
sions for the spherical functions in Eq. (C5) to obtain 



-In 



i^ip ■ 



f^ip - ko 



= -t^ipR- 



(C6) 



In this limit, the eigenvalues pip are therefore localized 
in the vicinity of a roughly circular line in the complex 
plañe given by 



^(l//i) 


- ko 


2 


^4ÍKÍl/fi)R 


tí{l/p) 


-f feo 







= 1. 



(C7) 



Let US now study the eigenvectors. Using standard 
properties of spherical harmonics and spherical Bessel 
functions [33 , we can show that 



- ji-i{i^ipR)ji+i{i^ipR) ^hi'^' 



ji{i^ipRf 



m^ — m' ^p^p' 



(C8) 



From the normalization condition {Cimp\Tli' m' p') — 
^i,i'^m,m'Sp^p', we find that A^p(r) = (-l)'^7^|(_^)p(r)* 
and 



Aip 



1 



Vji (f^ipR)^ - ji-Á^ipR)ji+i (i^ipR) 

On the other hand, we also have 



(C9) 



{"^Imp m' p' ) 



R^A^pAlp' 



^Ip' ^Ip 



f^ípjl-l{f^ípR)jl{f^lp^R) 



l^lp'jl-l{l^lp'R)jl{l^¡pR) 



(CIO) 



and {Cimp\¿^l'm'p') — {^lmp\^lmp')^l,U^m,m' ' It is UOW 

convenient to introduce a new coefficient 
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c, 



Ipp' 



K¡ Rjl-i{K¡ R)jl{Klp,R) - Klp,Rjl-i{Klp,R)jl{KrR) 



3i{K¡^Rf - ji.i{KlR)ji+i{^lR) 



Íl{Klp'RY - jl-lÍKlp'R)jl+lÍKlp'R) 



(Cll) 



in terms of which Eqs. (iCll) and (|C2|) become 



{21 + 1)11% 



9 ^~9^^lp 



\9? A^???(l 



[21 + l)nip'H¡pCipp' 



-f {1 - gij,ip,){l - gfiip)*' 



(C12) 



(C13) 



l P p' 

To find the borderline of the support of eigenvalue den- 



sity of the matrix (19) shown in Figs. [TJc) and (d), we 
apply the fohowing recipe. (1) Find solutions nip of Eq. 
(C5) numericaUy and then compute the corresponding 
¡iip. (2) Compute the coefñcients Cipp' using Eq. (Cll). 
(3) F ind hnes on the complex plañe g defined by Eq. 
( C13| ). (4) Transform the hnes on the complex p lañe g 
into contours on the complex plañe z using Eq. (C12). 



The latter contours are the borderlines of the support of 
eigenvalue density p(A). 
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